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C\J , Abstract 

o ■ 

psj ■ Random non-Hermitian Jacobi matrices Jn of increasing dimension n are con- 

sidered. We prove that the normahzed eigenvalue counting measure of Jn converges 
weakly to a limiting measure // as n ^ oo. We also extend to the non-Hermitian 
case the Thouless formula relating fi and the Lyapunov exponent of the second-order 
Q^, difference equation associated with the sequence J„,. The measure /i is shown to be 

^ I log- Holder continuous. 

> ■ 

S ■ 1 Introduction 

_> : 

^ ! Let ttj, bj, and Cj are three given sequences of complex numbers. Consider the second-order 
H ' difference equation for / 

ajfj-i + bjfj + Cjfj+i = zfj, j = 1,2,.... (1.1) 

This equation can be also written as 

jj.""' ]=gAy^),3 = 1,2,..., where g, = I 1 . (1.2) 



o 



Denote by fj{z) the solution of (jl.lj) satisfying the initial condition /o = 0, /i = 1. In 
terms of the transfer matrix Sn{z) = gn ■ . . . ■ gi, 

fn+l{z) \ _ c ,\ I 1 
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Obviously, fn+i{z) is a polynomial in z of degree n, 



n 



n 



fn+l{z) = knW{z - Zi), kn = Wl/Cj. 

1=1 j=l 



(1.4) 



Its roots the eigenvalues of the tridiagonal (Jacobi) matrix 



I hi ci 

02 &2 



\ 



n 



(1.5) 



a. 



-n— 1 
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In this paper we are concerned with the limiting distribution of the eigenvalues of J„ as 
ra — i> oo for random aj, bj, and Cj. 

If all bj are real and cj = a*_^i for all j the matrices J„ are Hermitian. The eigenvalue 
distribution of such matrices was studied extensively in the past in the context of the 
Anderson model, see e.g. PHIISI- In this case the eigenvalues are always real and there are 
several ways to prove that the normalized eigenvalue counting measure of Jn converges 
to a limiting measure as ^ oo. None of these proofs works in the non- Hermitian case 
and little is known about the limiting eigenvalue distribution of random non-Hermitian 
Jacobi matrices, however, see |H| lllj. 

Our interest to such matrices is partly motivated by non-Hermitian quantum mechan- 
ics of Hatano and Nelson CD] which, in one dimension, leads to equation ()1.1|) with the 
coefficients aj, bj, and Cj chosen randomly from the special class defined by the restrictions 



In this class the Liouville substitution^ reduces equation ()1.1|) to the symmetric equation 



tian case as the choice of boundary conditions to accompany equation (jl.lj) has a profound 
effect on the spectrum of the associated Jacobi matrix. If the Dirichlet boundary condi- 
tions, /o = and fn+i = 0, are chosen then the corresponding Jacobi matrix is Jn fll.^j) . 
As the Dirichlet boundary conditions are preserved by the Liouville transformation, the 
spectrum of J„ is real provided the coefficients {aj,bj,Cj) belong to the Hatano-Nelson 
class ()1.6|) . On the other hand, if one imposes the periodic boundary conditions, /o = /„ 
and /i = fn+i, then the spectrum of the corresponding Jacobi matrix turns out to be 
complex. This is not surprising of course as the Liouville substitution transforms the 
periodic boundary conditions for / into highly asymmetric boundary conditions for ip. 

1/^ = where 9i = 1 and 0k = (11^=! a*+i/cj)^^^ for k>2. 



bj G M and aj_^i/cj > for all j. 



(1.6) 
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What is surprising however is that in the hmit n ^ oo the complex eigenvalues lie on 
analytic curves J3] and are regularly spaced even if the coefficients in equation are 
chosen randomly ^Hj- These effects are specific to the Hatano- Nelson class and the proofs 
and analysis of the limiting eigenvalue distribution given in jT^l IE] exploit the relation 
between equations (jl.lj) and ()1.7j) . Of course, in the general case of arbitrary coefficients 
no such relation exists and one requires a different approach in order to investigate the 
eigenvalue distribution of J„. We develop such an approach in the present paper. 
Throughout this paper we assume that: 

Al {{cij,bj,Cj)}'^i is a sequence i.i.d. random vectors. 

A2 For some 6 >0 E[\aj\^ + \aj\-^ + \bj\^ + \cj\^ + \cj\'^] < oo. 

A3 The support of the probabihty distribution of the random vector (ai, &i, ci) contains 
at least two different points (a, b, c) and (a', b', c'). 

If all mass of the probability distribution of {aj,bj,Cj) is concentrated at one point 
(a, b, c) then of course we have a tridiagonal matrix with constant diagonals. This is 
a particular case of Toplitz matrices. Eigenvalue distribution of non-Hermitian Toplitz 
matrices was extensively studied in the past, see e.g. survey |2T] . 

Our main result expresses the limiting distribution of the eigenvalues of J„ in terms 
of the (upper) Lyapunov exponent 

7(^)= lim^log[|/„+i(z)p + |/„(z)p] 

of equation It is well known that (for every complex z) the above limit exists with 

probability one and is nonrandom. This follows from Oseledec's multiplicative ergodic 
theorem ^^I- A more subtle fact is that in our case 7(2;) can be calculated using the well 
known Furstenberg formula [Zj, and moreover 

7(^) = lim -Elog||^„(^)||, zee. (1.8) 

n— >oo n 

The function 7(2;) is subharmonic in the entire complex plane jl] and bounded from below, 

7(z) > ^Elog|ai/ci| for all (1.9) 

This inequality easily follows from det Sn{z) = YYj=i^j/(^j- The subharmonicity implies 
that A7, where A is the distributional Laplacian in variables Re 2 and Imz, defines a 
measure on C, see e.g. jT2]. Our main result is as follows. 

Theorem 1.1 Let fin be the normalized eigenvalue counting measure of Jn, i.e. Hn = 
^ '^^^i 5zi, where zi, . . . Zn are the eigenvalues of Jn - Then: 

(a) With probability one, fin converges weakly to fi = as n 00. 
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(b) (Thouless formula) For every z & C 

7(2)= / log \w — z\ dfi{w) — E log \ci\. (1-10) 
Jc 

(c) The limiting eigenvalue counting measure fi is log-Holder continuous. More pre- 
cisely, for any B^q^s = {z : \z — Zo\ < 6} , < 6 < 1, 

< (1.11) 

log I 

where C{zo, 6) ^ as 6 ^ 0. 



We deduce Theorem 11.11 from Theorem 11.21 which is of independent interest in the 
context of second order difference equations. 



Theorem 1.2 With probability one 

lim -log|/„+i(2;)| =7(2;) (1.12) 
for almost all z with respect to the Lebesgue measure on C. 



To prove Theorem II. 2[ we use the theory of products of random matrices. Of course 
this is unnecessary in the Hermitian case. In this case ()1.12|1 and the Thouless formula 
()1.1()|1 follow directly from the fact that fin converges weakly to a limiting measure fi 
m E] and the latter can be established independently and by more elementary means. We 
would like to emphasize that in the non-Hermitian case we follow the opposite direction 



route: the weak convergence of /i„ and the Thouless formula are deduced from ()1.12|) . To 
this end we make use of the relation 

/Wn = 77^Alog|/„+i|, (1.13) 
ZTrn 

where the equality is to be understood in the sense of distribution theory. Relation (jl.lHj) 
is well known in the function theory. It holds for arbitrary polynomial of degree n and can 
be easily derived with the help of the Gauss-Green formula. In this general setup it was 
shown by Widom j^Hl ETj that if the measures //„ for all n are supported inside a bounded 
region and in the limit n ^ 00 the function p„(z) = log \ z — w\dfin{w) converges to a 
limiting function p{z) almost everywhere in the complex plane then /i„ converges weakly 
io fi = 2^Ap. We shall need the following simple extension of this result to the case when 
the supports of fin are not necessarily bounded. 



4 



Let An be a (deterministic) sequence of square matrices of increasing dimension n, 
and 

Pn{z) = -log \det{An- Zln)\ = / log | W - 2 1 (w) , 

where is n x n identity matrix and yU„ = ^Apn is the normahzed eigenvalue counting 
measure of A„. Define 



n— >oo 



T/j = limsup / log \w\ dnn{ui), R>1. (1-14) 

\w\>R 



Proposition 1.3 Assume that there is a function p: C ^ [—00, +00) such that Pn{z) — * 
p{z) as n ^ 00 almost everywhere in C If ti < +00 then it follows that p is locally 
integrable, fj, = ^Ap is a unit mass measure, 

log \w\ dfi{w) < Ti < +00, (1-15) 

|t«|>i 

and the sequence of measures fin converges weakly to fi as n ^ 00. If, in addition, 
limij_»oo Ti?, = then we also have that 

p{z) = / log \w — z\dfi{w). (1-16) 
Jc 

Remark. In view of fll.lSp . the integral on the RHS in ()1.16p is a locally integrable 

function of z taking values in [—00, +00). 

For the sake of completeness, we give a proof of this Proposition in Appendix El 

In order to estimate the tails of eigenvalue distributions as required in the above 

Proposition II .HI we use the following inequalities^: 

n < limsup -^logdet(J„ + AnA*J, (1.17) 
and for any R> 1 and 6 > 

TR < — ^limsup-^trlogi+^(/„ + A„A:). (1.18) 

log rC n^oo Z ' n 

These inequalities can be derived with the help of Weyl's Majorant Theorem, for details 
of derivation see Appendix iBl 

^Note that logdet(/„ + AnA*„) = trlog(/„ + A„A;^). 
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Let us now return to the random Jacobi matrices J„. Straightforward but tedious 
calculations show'^ that 

1 " 

- tr logi+'(/„ + JnJ:) <^Y1 + /^l v^f)' ^here v,- = (a„ bj, 9), 

for some a,f3 > independent of Vj's and n. Therefore if the random sequence Vj is 
stationary and 

£;iog^+'^[l + |vi|2] < 00 for some (5 >0 (1.19) 

then the Ergodic Theorem asserts that with probability one the limits in ()1.17p (jl.l8|) 
are finite which implies ri < 00 and lim/j^oo tr = 0, as required in Proposition 11.31 
The assumptions of stationarity and ()1.19j) are less restrictive than assumptions A1-A3. 
However we are only able to prove Theorem 11.21 (which is the main ingredient to our proof 
of Theorem II. Ij) under these more restrictive assumptions. 



2 Products of random matrices 

Our proof of Theorem ll.2l makes use of several facts from the theory of products of random 
2x2 matrices. We list these facts below (Propositions l^?T] - 12. 3|1 . 

Let z/ be a probability distribution on the group G/(2,C) of invertible complex 2x2 
matrices and Qk be an infinite sequence of independent samples from this distribution. 

As before = ■ . . . • (?i for n = 1, 2, . . . . By P(C^) we denote the projective space 
on which every non-degenerate matrix g acts in a natural way. Let k be a probability 
measure on P(C^). We say that g preserves k if K{g^^.B) = n{B) for any Borel set B 
(here g.x is the result of the action oi g on x E P(C^)). By Gi, we denote the closure of 
the subgroup of Gl{2, C) generated by all matrices belonging to the support of u. We say 
that Gy preserves k if k is preserved by every g E Gy. 

Proposition 2.1 Let X^f^ > Ag"'' be the singular values of Sn- If 

Elog\\g\\ and Elog\det g\ are both finite (2.1) 
then with probability one the following limits 

lim -logAf) =7j, J = 1,2, (2.2) 

exist and are nonrandom. 



The limiting values 71 and 72 are called the Lyapunov exponents of the sequence S^- 



■^For any Hermitian matrix H = ||-ffjfe||"fe=i we have H < D = diag(di, . . . , d„) with dj = X!fc=i 
j = 1, . . . ,n. Therefore if / is a nondecreasing function then, by the Courant-Fisher minimax principle, 

trlog/(ff)<tr/(i?) = ELi/(rfj)- 
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Proposition 2.2 If in addition to condition \2.1]] . no measure n is preserved by Gy then 
the Lyapunov exponents of the sequence Sn are distinct, i.e. 71 > 72- 



Proposition 2.3 If condition \2.1\} is satisfied and no measure k is preserved by Gy then 

(i) For any unit vector x the probability is one that 

lim -log||S'„x|| = 7i. (2.3) 

n^oo Ti 

(ii) If in addition E{\\g\\^ + ||(^^^||^) < 00 for some 5 > then for any positive e there 
is a constant p{e) > such that uniformly in x, \\x\ \ = 1, 

Pro6 (|log I I - n7i I > en\) < e-"^(^\ (2.4) 

Remarks. 1. As all norms in are equivalent, the choice of norm in ()2.3j) and ()2.4j) 
is not important. However it is convenient to deal with the standard Euclidian norm. 

2. Propositions 12.11 - 12.31 are well known in the classical case of the real matrices, see 
e.g. [I71I2] for proofs of Propositions 12.11 and 12.31 and [TlIT^ for proofs of Proposition 12.21 
For complex matrices. Propositions 12.11 and 12.31 are proved in the same way as in fr7\ |2] . 
However, the proof of Proposition 12.21 is somewhat different from that given in [3 CH! • 
We shall now discuss the necessary changes which would allow the interested reader to 
reconstruct the proof in question simply by examining the one in ^^1- Namely, the main 
ingredient of this proof is the fact that the mapping g Tg, where 

(TJ) {x) = f{g-'x)\\g-'x\rf, 

defines a unitary representation of the group SL{m,M.) in Hilbert space L2{Sm,dl) with 
dl being the natural Lebesgue measure on the unit sphere Sm G M"\ (Obviously, we are 
interested in the case when m = 2.) 

In the case of the complex space the representation is defined by 

(TJ) {x) = f{g~'x)\\g-'x\\-"^, 

in Hilbert space L2{Sm, dl) with dl being again the natural Lebesgue measure on the unit 
sphere Sm G C". After that the proof proceeds in the way suggested in [T^ . 

3 Proofs of Theorems 11.11 and 11.21 

In order to be able to apply Propositions l2. 11 - 12. 3l we have to verify that under assumptions 
A1-A3 our matrices gj defined in ()1.2|) satisfy the conditions of these Propositions. 
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Is is apparent that assumption A2 guarantees that condition ()2.1|) is satisfied and 

remains to check that assumption A3 imphes that no measure 
K is preserved by G„ (here u is the measure induced on the group of matrices by the 
distribution of (ai, 61, ci)). To this end we note that if 

( z~b —a \ 

I Q J and ^' = 

then 

99' = rA" " . and ^' g=[ (..^v X ) and 



1 



a'c 



It remains to check that for almost all z that the group G generated by the matrices g, g' 
is rich enough in the sense that no measure is preserved by all matrices of this group. 
The main idea is as follows. For a "typical" z we construct two matrices, say B and D, 
from G such that the eigenvalues of B are of different moduli. It is easy to see then that 
the only measure preserved by all matrices of the form S", — 00 < n < 00 is the one 
supported by the lines in P(C^) generated by the eigenvectors of B. The matrix D & G 
is then chosen so that its action on P(C^) does not preserve these lines which means that 
the measure in question does not exist. We would like to emphasize that the presence of 
the parameter z plays a crucial role in this situation. 
More precisely, if z is such that 

2 arg(2; — h) ^ arg(ac) 

then the matrix g has eigenvalues with different moduli. In other words the moduli are 
different if z does not belong to a certain half line. The g' plays then the role of D (once 
again when z lies outside of certain curves). This statement can be checked by direct 
calculation and is sufficient for our purposes. 

However, in some important cases much more precise statements can be made. In 
particular if ^ = 1 then each of triangular matrices {gg'~^) and g'~^g is non-trivial for 
all but may be two values of z and a similar idea applies, see [2] page 213. 

Now we are in a position to apply Propositions 12. II - 12.31 For any two non-zero vectors 
X and y define 



d{x,y) 



{x,x){y,y)' 

where (■, ■) is the scalar product in C^. The function d{x, y) is the natural angular distance 
between x and y on the projective space P(C^). 

The following Lemma is the key element in the proof of Theorem II. 21 (In this Lemma 
and thereafter the abbreviation a.s. refers to the probability measure, i.e. any equality 
with the letters a.s. above it holds with probability one) 

Lemma 3.1 Suppose that the conditions of Propositions Wl\ - \2. 'A are satisfied. If yn is 
a sequence of random unit vectors in such that 

\\Snyn\ \ = e"'^^'*'^", where =' o(n) as n 00, (3.1) 
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then for any fixed unit vector x and any 5 > there is a constant r{x, 6) > such that 

Prob {d{x, Vn) < e-"-^} < e-""(^'^) (3.2) 

for all sufficiently large n. 

Proof. For any n one can always find two orthogonal unit vectors u„ and w„ such that 
SnSnUn = A^"^ and S^SnVn = ^2 ■ vicw of Propositiou \'2.1\ 

\\SnUn\ \ = e"^^+'" and \\SnVn\ \ = e"^2+^", where e'^,e'^ = o{n). 

Obviously the sequence Vn satisfies condition (jH.lj) and we first prove the large deviation 
estimate ()3.2|) for this sequence. 

Let a; be a fixed unit vector. Then x = {x, Un)un + {x, Vn)vn for every n, and, since 
= d{x,Vn) and |(x, f„)| < 1, we have that 

||S'„X|| < d{x,Vn)\\SnUn\\ + \\SnVn\\- 

Therefore if d{x, Vn) < e""'' then 

log < n7i + log(e-"^+^" + e-''^""-^^'>+''^) , 

and hence with probability one, 

log I |S'„x| I — nji < —n min((5, 71 — 72) + o(n). 

It follows now from Proposition 12.31 that 

Prob {dix, Vn) < e-"-^} < e-"'"^"'^) (3.3) 

for some r(x, 6) and all n > uq where hq depends on the matrices S'„, and also on x and 
5. 

Now, let yn be an arbitrary sequence of random unit vectors satisfying condition ()3.H) . 
and let y^ be a sequence of unit vectors orthogonal to i.e. {yn,yn) = for all n. 
Obviously, diun^y^) = K'^njl/n)! and, since S^SnUn = e'^^^^~^^''"Un, we have that with 
probability one 

d{un,y^) = e-2"^^+°(")|(5„u„,5„2/OI < e-"(^^-^^)+°("). 

It is then apparent that d{vn, yn) is also exponentially small for large n and therefore the 
large deviation estimate ()3.2|) for y„ follows from ()3.3|1 . □ 

Proof of Theorem \1.2l Let 

- = ™d """(/lit')' " = 
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Then 

.2. „ ^ ^ l/n+l(^)P ^ l/n+l(^)P 

"^^"^'^"^ |/«+lWP+|/nWP \\yn\? ' 

and therefore ^ ^ ^ 

-log |/n+l(2:)| = -logrf(x,y„) + -\og\\yn\\. (3.4) 

n n n 

In view of fjl.3p and Proposition I2.3n ). 

hm -log =• 7i(2r), (3.5) 

n— >cxD Tl 

where 7(2;) is the upper Lyapunov exponent of the sequence of transfer matrices Sn{z). 
On the other hand, S~^{z)yn = (1,0)-'", and therefore 

lim — log — = —71(2;). 

It follows now from Lemma l3.ll (applied to the matrices S~^{z) and the vectors x and 
yn/||yn||^) and the Borel-Cantelli Lemma that 

lim - log d{x, yn) =' 0. 

n— >oo n 

Therefore, in view of ()3.4|) and ()3.5p . for any fixed z the probability is one that 

lim - log \fn+i{z)\ = 71(2;). (3.6) 

n— >oo fl 

But then the probability is one that ()3.(ij) holds almost everywhere in the complex plane. 
This follows from the Fubini Theorem. Our proof of Theorem 11.21 is complete. 

Proof of Theorem li.il As explained in introduction, under assumptions Al - A3, 
the probability is one that Ti < C for some non-random C < +00 and lim^^oo tj? = 0. 
Therefore parts (a) and (b) of Theorem 11.11 follow immediately from Theorem 11.21 bv the 
way of Proposition 11.31 

The log-Holder continuity of /i is a corollary of the Thouless formula and the fact that 
the Lyapunov exponent 7(2) is bounded from below. This is very much in the same way 
as in the Hermitian case, see |3]. 

To prove (jLlip . we first note that the integral log \ w — z\dfi{w) converges absolutely 
for every z. Indeed, it follows from ()1.15|) that 

/ log \w — z\d^{w) < +00, 

J\w~z\>l 



*If 71 and 72 are the Lyapunov exponents of a sequence Sn then the sequence S'„ ^ has the Lyapunov 
exponents —72 and —71. 
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and this inequality together with the Thouless formula and the lower bound ()1.9p imply 
that 



as well. Therefore, 



/ \ log \w — z\\dfi{w) < +00 

J\w-z\<l 



C{z,S):= \ log \w - z\\dfi{w) ^ as(5^0. (3.7) 

J \w—z\<5 

Obviously, for 5 < 1, 

C{z, 6)= [ log(l/|ti; - z|)rf/iH > ^^^^''^ 



\<s " ^ - log(l/5)' 

and part (c) of Theorem II. II follows. Our proof of Theorem 11.11 is now complete. 

A Appendix 

Proof of Proposition \1. ,'A The local integrability of log|2;| and the condition ri < +00 
imply that the functions Pn{z) are uniformly integrable on bounded sets in C. It follows 
from this that p{z) is locally integrable and p„ — >• p as n — > 00 in -D'(C), the space of 
Schwartz distributions in C. Since A is continuous on distributions, we also have that 
Apn Aj9 in D'{C). Obviously Ap > 0, hence Ap is defined by a measure, see e.g. 

• As any sequence of measures converging as distributions must converge weakly we 
conclude that = ^Apn — *• /i = ^Ap weakly as measures. 

For any i? > 1, 

I djJniw) < ^ f log\w\dHn{w). 

J\w\>R ^Og\K\ J\u,\>l 

Therefore the inequality ri < +cxd implies that the sequence of measures /i„ is tight, and 
hence cannot lose mass. As each of fin has unit mass, so has the limiting measure /i. 
It follows from the weak convergence of /x„ to /i and ()1.14p that 



log < lim / loglwldfiniw) < ti 

w\<R J l<\w\<2R 

for any R> 1. This implies ()1.15p . Similarly, if lim ji^^tji = then 



lim / log \w\d jj,{w) = 0. (A.l) 



It remains to prove relation ()1.16|) . It will suffice to show that 



[ log\w - -Idiiiw) mD'{C) (A.2) 
Jc 
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when n — > oo. Let ip{z) be a continuous function with bounded support. Then 



c 



with 

g{w) = / 'iIj{z) log \w — z\(fz. 
Jc 

The function g is continuous and g{w) = 0(log|u;|) when \w\ — > oo. Assume now that 
hniij^oo tr = 0. Then 

hm hmsup / \g{w)\d^n{w) = 0, 

R^oo n^oo J\w\>R 



and 



hm / \g{w)\d^{w) = 

R—^oo 



'\w\>R 

because of flA.l|) . It now foUows from the weak convergence of to fi that 



hm / g{w)dnn{w) = / g{w)dn{w). 

Therefore 

hm / pn{z)'ijj{z)d''^z = I g{w)dfi{w) = I ■{ I log \w — z\dfi{z) \- ip{z)d'^z, 



c 



and ()A.2|1 foUows. 

B Appendix 

Derivation of inequalities and il.l^) . Let Zi, . . . ,Zn and Si, . . . , s„ be respectively 

the eigenvalues and singular values of An labeled so that \zi\ > \z2\ > ... > \zn\ and 
Si > S2 > ■ ■ ■ > Sn- Weyl's Majorant Theorem, see [Hj, page 39, asserts that 



3=1 i=i 



, ^, . . . , 77., 



for any function F{t) (0 < t < oo) such that F{e^) is convex on M. Obviously the function 
log'^''''^(t) satisfies this requirement for 5 > 0, and therefore 



m 

l+<5 



f log'^' \w\df,n{w) = -J2 log'^'k.l < - E log'^'l^.l < -Elog 
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where m is the number of eigenvalues of An such that \zj\ > 1. Obviously, 

m ^ m ^ n ^ 

j=l j=l j=l 

and therefore 

/ logi+^ \w\dfin{w) < -i^trlogi+^(/„ + AnAl), S > 0, 

J\w\>l ^ 

which implies fll.l7j) and ()1.18|1 . 
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